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^ ABSTRACT 
O 

^ Clouds of high infrared extinction are promising sites of massive star/cluster 

O formation. A large number of cloud cores discovered in recent years allows in- 

O vestigation of possible evolutionary sequence among cores in early phases. Wc 

have conducted a survey of deuterium fractionation toward 15 dense cores in var- 
ious evolutionary stages, from high-mass starless cores to ultracompact H ll re- 
gions, in the massive star-forming clouds of high extinction, G34.43+0.24, IRAS 
18151-1208, and IRAS 18223-1243, with the Submillimeter Telescope (SMT). 
Spectra of N2H+ (3-2), N2D+ (3-2), and C^^O (2 - 1) were observed to derive 
the deuterium fractionation of N2H+, Arac = A^(N2D+)/A^(N2H+), as well as the 
CO depletion factor for every selected core. Our results show a decreasing trend 
5^ in -Dfrac with both gas temperature and linewidth. Since colder and quiescent gas 

is likely to be associated with less evolved cores, larger Z^frac appears to correlate 
^ with early phases of core evolution. Such decreasing trend resembles the behavior 

psj of i^frac in the low-mass protostellar cores and is consistent with several earlier 

studies in high-mass protostellar cores. We also find a moderate increasing trend 
O of -Dfrac with the CO depletion factor, suggesting that sublimation of ice mantles 

O alters the competition in the chemical reactions and reduces Dfrac- Our findings 

^ suggest a general chemical behavior of deuterated species in both low- and high- 

^ mass proto-stellar candidates at early stages. In addition, upper limits to the 

ionization degree are estimated to be within 2 x 10^^ and 5 x 10~^. The four 
^ quiescent cores have marginal field-neutral coupling and perhaps favor turbulent 

cooling fiows. 

Subject headings: ISM: abundances — ISM: IRDC — ISM: individual (G34.43+0.24, 
IRAS 18151-1208, IRAS 18223-1243) — stars: formation 
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1. Introduction 

Clouds of high infrared (IR) extinction, such as infrared dark clouds (IRDCs), are dis- 
covered through silhouette against the bright, diffuse IR background emission of the Galactic 
plane (Egan et al. 1998; Rathborne, Jackson, & Simon 2006; Rygl ct al. 2010). Because of 
their cold (T < 20 K), dense {nu^ > 10^ cm^^), and massive (M > lO^M©) nature, IRDCs 
are thought to be promising sites of massive star/cluster formation. They often harbor cores 
in a sequence of evolutionary stages (Beuther & Sridharan 2007; Beuther et al. 2010), from 
fairly quiescent high-mass starless cores (HMSCs; Sridharan et al. 2005; Beuther & Srid- 
haran 2007), followed by cores with accreting low/intermediate- mass proto-stars/clusters 
(Beuther & Steinacker 2007; Wang et al. 2011), to high- mass protostellar objects (HMPOs; 
Sridharan et al. 2002; Beuther et al. 2002a), and even to ultra-compact Hll (UC Hll) 
regions (Battersby et al. 2010). In recent years, a large number of massive cloud cores 
(M > lO^M©) have been identified in IRDCs (Peretto & Puller 2010) and enable the in- 
vestigation of possible evolutionary sequence among cores in early phases (Chambers et al. 
2009; Battersby et al. 2010; Henning et al. 2010). In general, the strategy to probe evo- 
lutionary stages usually involves several indicators, each of which traces some star-forming 
activities, e.g. radio emission for ionizaiton, gas or dust temperature for internal heating, 
molecular linewidth for turbulence, asymmetry of line profiles for infall motions, or high- 
velocity gas for outfiows (e.g. Crapsi et al. 2005; Battersby et al. 2010). In our previous 
study of G28.34+0.06, we found a moderate decreasing trend in the deuterium fractionation 
of N2II+, -Dfrac = iV(N2D"'")/A^(N2H"'") , with evolutionary stage in three selected cores (Chen 
et al. 2010). Here we further investigate the possibility to use -Dfrac as an evolutionary probe 
for high-mass protostellar/cluster cores by searching its correlation with some indicators in 
three more clouds of high IR extinction. 

In early stages of star formation process, the low-temperature environment nurtures a 
peculiar chemistry with high deuterium enrichment as a result of exothermic deuteration 
reactions and significant depletion of CO and other neutral species. Deuterium enrichment 
is primarily initiated by the reaction 

H+ + HD ^ H2D+ + H2, (1) 

which is exothermic in the forward direction with AE/k = 230 K (Millar, Bennett, & Herbst 
1989). At very low temperature (T < 20 K), the back reaction becomes negligible, which 
results in an enhancement of 1120^^ (Stark, van der Tak, & van Dishoeck 1999; Stark et al. 
2004; CaseUi et al. 2003, 2008; Vastel et al. 2006; van der Tak, CaseUi, & CeccareUi 2005; 
Harju et al. 2006; Priesen et al. 2010) and even multiply deuterated H3 (Vastel, PhiUips, & 
Yoshida 2004; Parise et al. 2010; Roberts, Herbst, & Millar 2003). Deuteration is further en- 
hanced in cold cores after the removal of the gas-phase CO, the major destroyer of H2D+, due 
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to depletion of molecular species on dust grains (Ceccarelli et al. 2007; Bergin & Tafalla 2007; 
van der Tak 2006; Aikawa 2008; Millar 2005). Indeed, deuterated species are often observed 
with an enhancement of 2 — 3 orders of magnitude in star-forming cores (Crapsi et al. 2005; 
Fontani et al. 2006, 2011; Pillai et al. 2007) over the local interstellar value of 1.51 x 10~^ 
(Oliveira et al. 2003). Chemical models anticipate some correlation between the deuterium 
fractionation and the CO depletion factor, f^, which increases in the prestellar phase and 
declines later in the protostellar phase as the mantle sublimation occurs (Ceccarelli et al. 
2007; Crapsi et al. 2005; Emprechtinger et al. 2009). Such correlation has been recognized 
with the deuterium fractionation of N2H+ in a complied sample of low-mass prestellar and 
protostellar cores, among which cores in same molecular clouds show better correlation that 
indicates environmental influences, such as magnetic field strength, amount of turbulence, 
external radiation field, etc (Crapsi et al. 2005; Emprechtinger et al. 2009). On the other 
hand, early studies of high-mass protostellar objects often did not show a consistent behav- 
ior of deuterium fractionation with evolutionary stage other than a general enhancement 
(Fontani et al. 2006; Pillai et al. 2007). In the study of G28.34+0.06 (Chen et al. 2010), 
we started off our sample with cores in one IRDC and found a moderate decreasing trend in 
the deuterium fractionation of N2H"'" with evolutionary stages from a massive starless core 
(MM9), a more evolved core with fragmentation and outflows (MM4; Wang et al. 2011), 
to an UC Hll region (MMl). Subsequently, Fontani et al. (2011) significantly improved 
the statistics with a larger sample of twenty-seven cores and also revealed a similar trend 
across evolutionary phases from HMSCs, HMPOs, to UC Hll regions. In particular, the 
authors reported a nice decreasing trend with gas temperature and column density of N2H^. 
Recently, Miettinen, Hennemann, & Linz (2011) reported a similar decreasing trend but in 
the A^(DCO^)/A^(HCO^) ratio with gas temperature. They further made first estimates of 
ionization fraction and cosmic ionization rate in massive IRDCs. Over all, the use of the 
A^(N2D+)/A^(N2lI+) ratio as an evolutionary probe to high-mass protostellar candidates has 
been suggested (Chen et al. 2010; Fontani et al. 2011; Caselh 2011; Miettinen et al. 2011). 

We select cores in three massive star-forming regions of similar luminosity (L ~ 10^ Lq) 
to probe a sequence of evolutionary stages: one IRDC, G34. 43-1-0. 24 (hereafter G34.43; Rath- 
borne et al. 2006), and two HMPOs, IRAS 18151-1208 and IRAS 18223-1243 (hereafter 
118151 and 118223, respectively; Sridharan et al. 2002). All three regions harbor multiple 
cores at different evolutionary stages, and discussions of individual clouds are provided in 
the Appendix. To better parameterize the evolutionary stage with gas properties, we limit 
our sample to cores with available ammonia gas temperatures from a survey conducted by 
Sakai et al. (2008). This renders eight cores in 034.43, three cores in 118151, and four 
cores m 118223. Table [l] lists the selected cores and their known properties such as gas 
temperature, Tg, integrated flux density at 1.2 mm, Fi.2mm, core size, 2R, molecular number 
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density, riHa- In addition, four out of the fifteen selected cores are not associated any known 
outflows (Table [T| and possibly at the earliest stage. We consider them to be quiescent and 
perhaps pre-stellar cores in this study. Our sample of cores have gas temperatures from 14 
to 21 K, which are within a promising range to detect variation in /£> based on dust mantle 
sublimation (CoUings et al. 2003). These cores are generally colder than those considered 
in previous studies (Fontani et al. 2006, 2011). 

Based on previous studies, e.g. Chen et al. (2010) and Fontani et al. (2011), we 
would like to assess the use of -Dfrac as an evolutionary probe with a large sample of cores, 
particularly for evolutionary stages before HMCs and UC Hll regions. In this study, a 
similar strategy is employed to compare -Dfrac among cores in similar IRDCs to reduce the 
environmental fluctuations among selected cores. Moreover, we also investigate the behavior 
of the A^(N2D+)/A^(N2H"'") ratio with the CO depletion factor through C^^O observations 
and estimate the electron abundance. In § [2| we describe the observations and related 
parameters. We then explain in § [3] the derivation of Dfrac with self-consistent spectral fits 
for N2H"'" and N2D"'" as well as the determination of the CO depletion factor. Lastly, in § |4| 
we discuss the use of the iV(N2H"'")/A^(N2D"'") ratio as an evolutionary probe as well as the 
upper limits for the degree of ionization and its implications. 



2. Observations and Data Reduction 

We observed N2H+ (3 - 2) at 279.511780 GHz, N2D+ (3 - 2) at 231.321864 GHz, and 
C^^O (2 - 1) at 219.56036 GHz towards selected cores in G34.43, 118151, and 118223 with 
the 10-meter Arizona Radio Observatory (ARO) Submillimeter Telescope (SMT) on Mount 
Graham, Arizona. The observations were carried out between November 2008 and April 
2011 in the beam-switching (BS) mode for N2H"'" and N2D+ and in the absolute position- 
switching (APS) mode for C^®0 to correctly subtract the background because C^^O emission 
is supposed to be extended. The pointing centers are given in Table [l] with typical uncer- 
tainty of ~ 3". The primary beam is about 27" for N2H"'", 32" for N2D"'", and 34" for 
C^^O. The spectral resolution is 1 MHz, corresponding to a velocity resolution of 1.07 and 
1.30 kms~^ for N2H"'' and N2D"'", respectively, and 0.25 MHz, equivalent to 0.34 kms~^, 
for C^^O. The temperature scale was obtained using standard vane calibration, and 
the main beam temperature, T^b, was derived through Tmb = T^/Vmh with a main beam 
efficiency r^mb = 0.75. Typical system temperature during the observations was around 
200 — 300 K and the respective rms noise level for each observation is given in Table [Tj 
Data reduction was performed with the CLASS package (Guilloteau & Lucas 2000, ; see also 
http://www.iram.fr/IRAMFR/GILDASl. 
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3. Data Analysis and Spectral Fits 

Except the N2D"'' in G34.43-MM8 (Fig. [l^), all three molecular lines are detected in 
emission toward every core. The spectra for cores in G34.43 are shown in Fig. [T] while those 
for in 118151 and 118223 in Fig. [2]and[3| respectively. We derive -Dfrac by analyzing the N2H"'" 
and N2D"'" spectra with a self-consistent spectral fit. The CO depletion factor, fo, is also 
estimated with the observed C^^O spectra and available 1.2 mm continuum flux density. 



3.1. Self-Consistent Spectral Fits for N2H+ and N2D+ 

In some cores, the emission of N2D+ is undesirably weak and becomes more challenging 
to constrain the spectral fit parameters. Instead of performing an independent spectral 
fit for N2D+, we improve the model described in Chen et al. (2010) to perform a self- 
consistent fit for the N2H"'" and N2D+ spectra using -Dfrac a scaling factor. Because of 
the numerous hyperfine components, we fit each spectrum of every core with a synthetic 
spectrum comprised of thirty-eight hyperfine components with updated line frequencies and 
spontaneous emission rates (Pagani, Daniel, & Dubernet 2009). For each individual source, 
all the hyperfine components of every J-level are assumed to be in thermal equilibrium at 
a single excitation temperature, Tg, adopted from the ammonia gas temperature in Sakai 
et al. (2008). The synthetic spectra are described by three more parameters: the total 
column density, A^(N2H"'"), the systemic velocity, i^lsr? and the full-width at half maximum 
(FWHM) as line width. At;. Model spectra are optimized with the minimization of the 
reduced value, x^, and the results are listed in Table [2] Note that the uncertainties of the 
adopted Tg (Sakai et al. 2008) are also incorporated into the uncertainties of all the derived 
spectral parameters. 

Taking the cosmic background temperature, Tbg = 2.7 K into account, we derive the 
optical depth of an observed spectrum with 



t{v) = — In 



J{Tg)-J{n,)\ 



(2) 



where Tmhiv) is the main beam temperature of the spectra, and J{T) = {hi'/k)/{e^^/^^ — 1). 
In general, the spectra have fairly small optical depths if a beam filling factor of unity 
is assumed. The optimized spectral model also delivers an estimate for the N2H+ optical 
depth by integrating optical depths of all the hyperfine components. The maximum optical 
depth, Tmax, is about 0.8 in G34.43-MM1 (Table |2|). Emissions of N2H+ and N2D+ in all the 
observed cores remain optically thin. Although sub-structures within our observing beam 
sizes cannot be completely ruled out, the emission in these early phases of core evolution 
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is mostly attributed to large scales. In the example of G28.34+0.06, the integrated flux 
observed with the Submillimeter Array (SMA) is less than 10% of the total flux observed 
with the single-dish telescope SMT (Chen et al. 2010). 



3.2. The CO Depletion Factor 

The CO depletion factor, fo, is defined as the ratio of the canonical CO abundance, 
x(C0)can5 to the observed CO abundance, x(C0)obs5 

_ a;(CO)can _ 3;(C^®0)can /„n 
" x(CO)obs " x(C180)obs' ^ ' 

Using the C^^O abundance of 1.7 x 10~^ in the solar neighborhood (Frerking, Langer, & 
Wilson 1982) and the abundance gradients of A log[C/H]/Ai? = —0.066 dexkpc"^ and 
Alog[0/H]/Ai? = -0.065 dexkpc'^ in the Galactic Disk (Wilson & Matteucci 1992), the 
canonical abundance of C^^O is estimated to be 

a;(Ci«0)ean = 1-7 X 10"^ ■ lO-^'^^^ (^^GC-I^e)^ (4) 

where -Dqc is the Galactocentric distance of the core, and Dq = 8.5 kpc is the distance of 
the Sun to the Galactic Center. Given the location of our IRDCs, x(C^^O)can = 3.80 x 10^^, 
3.92 X 10-^ and 4.68 x 10"^ for G34.43, 118151, and 118223, respectively. 

Over all, the observed brightness temperature, Tmb, is smaller than the kinetic temper- 
ature, Tg, and renders fairly small optical depths. The maximum optical depth estimated 
with Eq. in individual IRDC is 0.5 for G34.43-MM2, 0.5 for I18151-MM3, and 0.8 for 
I18223-MM1. Assuming that all rotational levels are thermalized, we determine the column 
density of C^^O with the method based on Caselli et al. (2002) that accommodates the effect 
of background emission at Tbg, 

'^(^"O' - l!^^4; "'(^''O) S i^S^y 

where Q{Tg) is the partition function, /i^ 5*21 = 0.02440 Debye^ for the J = 2 — 1 transition, 
uq the transition rest frequency, E^p = 15.8 K the upper level energy, and W{C^^O) the 
integrated brightness temperature in velocity (Table [T|. Since some of the spectra do not 
resemble a Gaussian profile, direct integration in the channels with significant emission is 
performed to obtain W{C^^O) without fitting a Gaussian profile, and the derived A^(C^^O) 
is listed in Table [3l 
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The observed C^^O fractional abundance, x(C^^O)obs; depends on the column density 
of H2, and 2;(C^^0)obs = ^(C"'^^0)obs/^H2- In attempt to reduce the uncertainty 
in deriving x(C^®0)obs, we first match the angular resolutions between the C^'^O and dust 
continuum observations by convolving the 1.2 mm continuum maps in the literature (beam 
size = 11"; Beuther et al. 2002a; Rathborne et al. 2006) with the 34" beam of our C^^O 
observations. The column density of H2 is estimated with the 1.2 mm peak flux density, 
>S'i.2mm, arising from warm dust 

c 1 

AT ^1.2mm /„N 

where Ki.2mm = 0.005 cm^g~^ is the dust opacity assuming a gas-to-dust mass ratio of 100 
(Shepherd & Watson 2002), By{Td) the Planck function at dust temperature, T^, VLb the solid 
angle subtended by the convolved beam size of 34", = 1.36 the mean molecular weight, 
and the mass of H2 molecule. 

For warmer cores in our sample, the averaged dust temperature, T^, have been found 
from studies of the spectral energy distributions (SEDs) (Table 3; Rathborne et al. 2005; 
Marseille et al. 2008; Beuther et al. 2010). However, most cores in our sample remain 
as extinction features in the near- or mid-IR, and it is challenging to derive their dust 
temperatures. Sensitive mid- and far-IR observations, such as Herschel, will offer better 
opportunities to constrain dust emission properties, including temperature, in cold clumps of 
IRDCs (e.g. Beuther et al. 2010; Stutz et al. 2010). Alternatively, we adopt gas temperature, 
Tg, for the conversion between S'i.2mm to A^H2 when dust temperature is unavailable. In case 
the thermal coupling between dust and gas is not ideal, this assumption may underestimate 
Td and hence overestimate iV(H2) by a factor of < 2 in our sample. In the case of G34.43- 
MM2, no SED study is found to constrain its dust temperature. Since the core is associated 
with an UC H ll region, we expect its averaged dust temperature to be warmer with respect 
to MMl and sets a lower limit of > 34 K for MM2. 

Once the observed C^^O abundance, x(C^^O)obs, is determined, the CO depletion factor, 
fo, can be computed with Eq. ^ accordingly, and the results are summarized in Table |3} 



4. Results and Discussion 

Except G34.43-MM8 with undetected N2D"'" emission, all our cores show a general 
enhancement of 2—3 orders of magnitudes in -Dfrac (Table [2]) over the local interstellar value 
of 1.51 X 10~^ (Oliveira et al. 2003). Note that the -Dfrac in G34.43 spans a fairly large range 
from 0.0039 to 0.11, nearly a factor of 30. The deuterium fractionation, -Dfrac, are compared 
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with the gas temperature, Tg, the fitted hnewidth, At>, and the CO depletion factor, 
(Fig. |4|. A clear decreasing trend in -Dfrac with both Tg and Av but a weaker increasing 
trend with fu can be seen. Although these behaviors generally agree with expectations based 
on chemical models, an analytical formula to describe the dependence is not obvious with the 
large scatters in Fig. |4j To search for dependence, we perform statistical tests between Z^frac 
and other parameters to evaluate the correlation along with the significance, which gives the 
likelihood for the correlation occurring by chance. The decreasing trend between -Dfrac and 
Tg (Fig. [4^) has Spearman's p rank correlation coefficient of p = —0.67 with significance of 
0.6% and Kendall's r rank correlation coefficient of r = —0.50 with significance of 1.0% 
while the correlation between -Dfrac and Av (Fig. ^p) gives p = —0.61 with significance of 
1.6 % and r = —0.49 with significance of 1.2 %. On the other hand, the dependence between 
Dfrac and fi) shows larger scatters (Fig. [Ij:). When excluding G34.43-MM1 with an unusually 
large fo, we find an improved correlation that has p = 0.49 with significance of 7.8% and 
r = 0.34 with significance of 9.0 %. The concerns to include the HMC G34.43-MM1 will be 
elaborated in Sect. 14.21 

4.1. Deuterium Fractionation as an Evolutionary Probe 

Over all, a monotonically decreasing trend in -Dfrac with both increasing gas temperature, 
Tg, and fitted linewidth, Av, is discerned (Fig. |4^ and |4jD). While examining cores in 
individual clouds, there seems to have a shghtly better correlation, implying possible cloud- 
to-cloud variation due to the infiuence of environments as previously suggested in the studies 
of low-mass cores (Crapsi et al. 2005; Emprechtinger et al. 2009). Miettinen et al. (2011) 
also suggest relatively large environmental variations in the cosmic-ray ionization rates in 
massive IRDC cores. In particular, quiescent cores with no outfiow activities, i.e. G34.43- 
MM6, -MM9, I18151-MM3, and I18223-MM2, all have the lowest temperature (Table Q 
and the largest D^ac (Table |2| as well. Since warm and turbulent gas is more likely to 
be associated with evolved cores, the observed D^ac suggests a decreasing dependence with 
evolutionary stage. For a better determination of their evolutionary stage, one may desire to 
compare with theoretical evolutionary models, which depend on several physical parameters, 
such as bolometric temperature, bolometric luminosity, and envelope mass (Froebrich 2005). 
However, most of our selected cores, especially those with the largest Dfrac, do not have 
observations in the mid- and far-IR wavebands to better constrain their dust temperature and 
bolometric luminosity. Alternatively, we compare Dfrac with gas temperature and linewidth 
instead of bolometric temperature and luminosity. 

A few previous studies have reported an anti-correlation between deuterium fraction- 
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ation and evolutionary stage of massive proto-stellar cores. Chen et al. (2010) found 
-Dfrac = 0.017 — 0.052 in three cores at different evolutionary stages within the IRDC 
G28.34+0.06, including a massive starless core (MM9), a core with fragmentation and out- 
flow activities (MM4; Wang et al. 2011), and an UC Hll region (MMl), and suggested 
a decreasing trend in -Dfrac with evolutionary stage. A subsequent study by Fontani et al. 
(2011) significantly improved the statistics with a sample of twenty-seven cores and also re- 
vealed this decreasing trend with values of -Dfrac in the range of 0.012 — 0.7, 0.017— < 0.4, and 
0.017 — 0.08 for their observed HMSCs, HMPOs, and UC Hll regions, respectively. In par- 
ticular, they found an anticorrelation between -Dfrac and Tg. Miettinen et al. (2011) further 
reported a decreasing trend in the ratio of A^(DCO+)/A^(HCO+) = 0.0002 - 0.014 with gas 
temperature in their sample of seven IRDC cores. They also reported -Dfrac = 0.002 — 0.028 
for the four cores with higher gas temperature. Our values of -Dfrac are in the range of 
0.0039 — 0.11, which are comparable to the values obtained by Fontani et al. (2011) and 
Miettinen et al. (2011). The anti-correlation between -Dfrac and Tg is also seen in our results 
(Fig.§i). 

Given the gas temperature range of Tg = 14 — 21 K, the corresponding thermal linewidth 
of N2H"'" is merely Afth = 0.15 — 0.18 kms~^. The observed linewidth in the range of 2.3 — 
4.8 kms~^ is dominated by non-thermal motions, possibly arising from turbulent motions 
among clumps. In Fig. [sj we compare the linewidth of our N2II+ (3 — 2) spectra with 
linewidths of N2H+ (1 — 0) and NH3 (1, 1), (2, 2), and (3, 3) spectra observed by Sakai et al. 
(2008). The N2H+ (1 ~0) observations had a smaller beam of 18" while the NH3 observations 
had a much larger beam of 73" with respect to our 27" beam for N2H'*" (3 — 2). Between 
the two transitions of N2H"'", the J = 3 — 2 transition with higher -Eup = 26.8 K shows 
broader linewidth than the J = 1 — with lower E^p = 4.5 K as expected for warmer gas 
to be more turbulent (Fig. |5^). In general, NH3 lines have a much lower critical density of 
rzcrit ~ 2 X 10^ cm~^ (Evans 1999) and tend to trace the outer part of the cores with respect 
to N2H"'" lines with ncrit ~ 10^ cm~^ . Except in the two quiescent cores, G34.43-MM6 
and I18151-MM3, linewidths of the NH3 (1, 1) and (2, 2) spectra are comparable to those of 
N2H+ (1 — 0) but smaller than those of N2H+ (3 — 2) (Fig. [sjo andjsj^). The two quiescent 
cores are probably in a very early stage where turbulence dissipation may occur to produce 
smaller linewidth in the inner part of higher density (Goodman et al. 1998). In a number of 
more evolved cores, the NH3 (3, 3) emissions with much higher -Eup = 124.5 K are detected 
and show much larger linewidth, even up to 7.2 kms~^ in the case of G34.43-MM3 (Fig. |5]i). 
Since all these evolved cores are associated with outflow activities, the NH3 (3, 3) emission 
is tracing the hot gas which could be in the outflows (Zhang et al. 2002). 

Unlike NH3 and NH2D, which are affected by grain surface reactions (Giirtler et al. 
2002; Bottinelli et al. 2010), N2D"'' and N2H"'" are pure gas-phase reactants and do not 
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participate condensation and subsequent sublimation of ice mantles. Compared to other 
molecular species, the deuterium fractionation of N2H"*" better reflects the physical conditions 
at present time without being confused by evaporation of mantles which had formed at earlier 
times with enhanced deuteration (Emprechtinger et al. 2009). In a sample of Taurus cores, 
a clear increasing trend in the deuterium fractionations of both NH3 and N2H"'" was observed 
in prestellar cores (Crapsi et al. 2005; Hatchell 2003), whereas in protostellar cores, -Dfrac 
shows a faster decreasing trend with dust temperature than does the A^(NH2D)/iV(NH3) ratio 
(Emprechtinger et al. 2009). For high-mass protostellar cores within similar gas temperature 
range, we note that the correlation between -Dfrac and Tg appears stronger in our cores than 
does the A^(NH2D)/A^(NH3) ratio (Pillai et al. 2007). 



4.2. Deuterium Fractionation and the CO Depletion Factor 

To further examine the relationship between deuterium fractionation and the CO deple- 
tion, we computed the C^^O column density and the CO depletion factor, fc, as described 
in Sect. 3^ The results are listed in Table [3] and shown in Fig. |4|:. When excluding G34.43- 



MMl, which has the largest value of fo, we find an increasing trend in -Dfrac with f^- This 
agrees with the general expectations from chemical models that CO is the major destroyer 
of H+, H2D+, N2H+, and N2D+ (CaseUi et al. 1998). As the envelope heats up, the CO 
abundance is expected to quickly rise up based on the dramatical drop of the CO sublima- 
tion timescale from 10^ yr at ~ 12 K to 0.1 yr at ~ 20 K (Colhngs et al. 2003). The 
warmer temperature together with the return of gas-phase CO can alter the competition in 
the chemical networks and brings a drop in -D^ac (Roueff et al. 2005; Aikawa et al. 2005; 
Aikawa 2008). When CO returns to the gas phase, it will quickly react with N2H~*' and N2D"'" 
to form HCO+ and DCO+, respectively (Lee, Bergin, & Evans 2004). However, one should 
be cautious about possible chemical stratification once the CO sublimation starts in the cen- 
tral warm region. As a core warms up, the N2H"'" and N2D+ emissions tend to trace the cold 
outer region whereas the CO emission is dominated in the central region. In the study of 
the Ophiuchus B2 core, Friesen et al. (2010) found that -D^ac increases at greater projected 
distances from the embedded sources. When the observed emission arises from a partially 
filled volume, the beam-averaged abundance for each species will start to deviate from the 
actual abundance used in chemical models. Observations with high angular resolutions will 
help to image the spatial distributions and reduce the confusion in comparing abundances. 

In the case of G34.43-MM1, an unusually large value of A^H2 and hence fo is obtained. 
Millimeter interferometric studies (Cortes et al. 2008; Rathborne et al. 2008) reveal a very 
strong (-Lboi ~ 2 X 10^ L0) and compact (2-R ~ 0.03 pc) source that exhibits signatures of 
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an HMC, which is thought to have a typical temperature closer to 100 K. This source has 
developed a steep temperature gradient, from 34 to 100 K across core scales from 0.1 to 
0.015 pc (Rathborne et al. 2008), translating to a single power-law dependence of r"'^'^''. 
This steep temperature gradient suggests the presence of an inner region where optical depth 
is large to the IR photons that the photon diffusion should be considered (Kenyon, Calvet, 
& Hartmann 1993; Osorio, Lizano, & D'Alessio 1999; Chen et al. 2006). In our approach 
to estimate /d, the dust emission is likely dominated by the hot inner region while the CO 
emission arises from a cold and large envelope. The dust temperature, T^, derived from 
the SED fit depends on the fiux densities in the mid-IR wavebands that may suffer from 
significant optical depth and does not refiect the physical conditions of the hot central region 
where the peak of the optically thin 1.2 mm emission is produced. Such steep temperature 
gradient may cause an underestimate in and overestimates in N-^^ and fn in our current 
calculation. To avoid potentially misleading interpretation, we exclude the result of G34.43- 
MMl when discussing the electron abundance in the following section. 



4.3. The Ionization Degree 

The enrichment of primary deuterated ions, e.g. H2D"'", CH2D"'", and C2HD"'", will give 
rise to the enrichment of subsequent deuterated species, such as N2D+, but with lower [D]/[H] 
abundance ratios due to the statistical nature of the fractionation process. In simple steady- 
state models based on gas-phase ion-molecular chemistry, an upper limit to the electron 
abundance, be found with assumptions that all the deuterium enrichment originates 

in II2D+ and that the recombination on negatively charged grains is negligible (Wootten, 
Snell, & Glassgold 1979; Caselli et al. 1998; Caselli 2002). Following the method described 
in Caselli (2002), the deuterium fractionation, -Dfrac, may be expressed as a function of Xe 
and abundances of HD, a;(HD), and important neutral species, x{m), 

A,.^ (7) 

m 

where /chd = 1-5 x 10~^ cm^s~^ is the rate coefficient for the reaction in Eq. Q, = 
6 X 10^^ (T/300)^'^'^^ cm^ s~^ the dissociated recombination rate of H2D"'" (Caselli et al. 
1998), km the destruction rate for II2D+ due to reactions with neutral species m, such as 
CO and O. The numerical factor of 1/3 accounts for the statistical branching ratio of 1/3 
to transfer the deuteron in the reaction of II2D+ with N2. The HD abundance is taken from 
the interstellar value of x(HD) = 2[D]/[H] = 3 x 10"^ (Oliveira et al. 2003; CaseUi et al. 
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1998). The electron abundance is then given by 



^HD a;HD 1 

X. = — — — > Krry.Xim) 



3 ke -Dfrac 

m 

Since CO is the dominant neutral species that destroys 1320"^, we make an approximation 
by neglecting other ion-neutral reactions to get 



^kmx{m) > /i;coa;(CO) = kco 



x{CO), 



D 



(9) 



where kco = 6 x 10"^*^ (T/SOO)"*^'^ cm^s~^ is the H2D+ destruction rate due to reactions 
with CO (Caselli et al. 1998), and a;(CO)can = 1-5 x 10""^ is the canonical CO abundance at 
the locations of our cores. We may set an upper limit for x^ to be 

^HD a^HD 1 ^CO2;(C0)can 1 /-iriN 

~ ~^u~T) k T 

O Ate -f^frac l^e J D 



3.8x10'^ 9.7x10"^ 



-Dfrac / 



D 



(for Tg = 16 K). (11) 



For easy comprehension of the dependence on -Dfrac and fn, the numerical values are provided 
for the case of Tg = 16 K. With the derived -Dfrac and Jd, we obtain the ionization degree in 
the range of Xe = 2 x 10~^ — 5 x 10~^ for our selected cores (Table [s]). These values lie at high 
end of the ionization degrees reported in early studies of low-mass dense cores (Caselli et al. 
1998; Williams et al. 1998) and massive cores (Bergin et al. 1999). Recently, Miettinen 
et al. (2011) derived first estimates for ionization degrees in IRDC cores with upper limits 
of Xe = 2 X 10~^ — 2.9 X 10"'^ and lower limits of = 3 x 10~^ — 5.6 x 10~^. Our estimates 
give smaller values compared to their upper limits but are within the range bracketed by 
their upper and lower limits. The smaller Xe upper limits are mainly attributed to the larger 
deuterium fractionation derived from N2H+ instead of HCO"*". Furthermore, our Xe values 
show a moderate correlation with the evolutionary stage. Since more evolved cores show 
smaller values of I^frac and fn, the plausible correlation between Xe and evolutionary stage 
may lead to the decreasing trend of -Dfrac- Most of the cores in our sample have shown 
outflow signatures and likely have begun to form clusters with low- and intermediate-mass 
protostars. This increasing degree of ionization could be arising from accretion shocks and/or 
heating of the gas related to the central proto-stellar /cluster objects (Stabler, Shu, & Taam 
1980; Hosokawa, Yorke, & Omukai 2010; Calvet et al. 2004). Given the porous nature of 
the surrounding medium (Indebetouw et al. 2006), it is possible for part of the energetic 
photons to reach outer part of the core and rise up the overall electron abundance (Kim & 
Koo 2001). 
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The four quiescent cores, G34.43-MM6, -MM9, I18151-MM3, I18223-MM2, also have 
the lowest ionization degree, Xf. — 2x 10~^ — 9 x 10~^, among cores in same IRDC. In these 
cores, the degree of ionization may play a role in regulating star formation efficiency through 
ambipolar diffusion, in which magnetic fields drift relative to a background of neutrals (Mestel 
& Spitzer 1956; Shu, Adams, & Lizano 1987). In a partially ionized medium, charged 
particles are coupled to magnetic fields while neutrals are supported against their self-gravity 
through the frictional drag that they experience when drifting through ions. In addition to 
ambipolar diffusion, evolution of massive star-forming cores can be affected by various effects 
such as turbulence (McKee & Tan 2003), rotation, and magnetic fields. For a more complete 
picture, one needs to consider all the important supports against gravity (Myers & Goodman 
1988; McKee & Ostriker 2007). Here we just make an attempt to estimate the characteristic 
scale for ambipolar diffusion, ^adj which gives the smallest scale for which the magnetic 
field is well coupled to the bulk of the gas. Following the picture described in McKee & 
Ostriker (2007) and approximating the ion abundance with rij = XgnHa, we estimate the 
characteristic ambipolar diffusion scale to be 

Va r. r.. ( Va \ ( Ui 



f,,= ^ « 0.06 PC (^^)( 

rij aj„ \6 kms~^/ V 



10~^ cm~^ 



where va = B / ^/47fJ^m}l~n^ is the Alfven speed associated with the large-scale magnetic 
field B and aj„ ^ 2x 10^^ cm'^ s^^ the ion-neutral collision rate coefficient (Drainc, Robcrge, 
& Dalgarno 1983). For the unknown magnetic field strength in supersonic cores, we take the 
median value given by McKee & Ostriker (2007) (originally from Crutcher (1999)), 



where Av = y' Av'^ — Av^^ is the line width for nonthermal motions. Applying the values 
for and Av in the four quiescent cores, we find the magnetic field strength to be in a 
range from 210 to 572 fiG. The characteristic scale for ambipolar diffusion is in a range from 
£ad = 6 X 10~^ to 4 X 10~^ pc, much smaller than the typical size of the cores. 

Since most star-forming cores have Alfven speeds comparable to the nonthermal com- 
ponent of the velocity dispersions, the turbulence in these cores should have a substantial 
magnetohydro-dynamic (MHD) component (Myers & Goodman 1988; Crutcher et al. 1999). 
Myers & Lazarian (1998) have proposed a pressure-driven cooling flows associated with local 
dissipation of turbulence due to wave damping by ion-neutral friction in an inner core re- 
gion. From a sufficiently turbulent outer region, the MHD wave power transmission, g, into 
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a spherical inner core is described in terms of the field-neutral coupling parameter, W, which 
is defined as the ratio of the core size, R, to the minimum cutoff wavelength, Aq = tt^ad; for 
the propagation of MHD waves (Myers & Khersonsky 1995). Using the upper limits of Xe 
in the four quiescent cores, we obtain 14^ = 23 — 119 with a mean value of W = 75, roughly 
in the regime of marginal field-neutral coupling {g ~ 0.4 for W < 100; Myers & Lazarian 
1998). The MHD waves excited in the outer region would have a fairly limited range of 
allowed wavelengths to transmit the turbulence power into the dissipating inner regions, and 
a pressure-driven turbulent cooling fiow may occur in these cores. For comparison, Bergin 
et al. (1999) reported W — 20 m massive cores with ionization degree inferred from chemical 
models while Miettinen et al. (2011) obtained W < 18.5 using their lower limits of ionization 
degree. But one should be cautious about a direct comparison among these values because 
of different methods to estimate the ionization degrees and the magnetic field strength. 

5. Summary 

We observed emissions of N2H+ (3-2), N2D+ (3-2), and C^^O (2 - 1) toward 15 cores 
in the IRDC G34.43, and the HMPO 118151 and 118223. The main findings are summarized 
as follows: 

1. A clear decreasing trend in the deuterium fractionation of N2H+, -Dfrac, with evolution- 
ary stage traced by increasing gas temperature and hnewidth. This decreasing trend 
agrees with the findings in previous studies by Chen et al. (2010), Fontani et al. (2011), 
and Miettinen et al. (2011). An increasing trend, though with larger scatters, in -Dfrac 
with the CO depletion factor, /d, is also found. Such trend resembles the behavior of 
-Dfrac in the low-mass protostellar cores and suggests the use of the A^(N2D+)/A'"(N2H+) 
ratio as an evolutionary probe to high-mass proto-stellar/cluster candidates. 

2. A significant enhancement of 2 — 3 orders of magnitude in the A'"(N2D+)/A'^(N2H+) 
ratio in all detected sources over the local interstellar [D]/[H] ratio. Such enhancement 
agree well with those observed in other massive star-forming cores. 

3. The upper hmits of electron abundance are estimated to be in the range from 2 x 10~^ 
to 5 X 10~^, which lie at the high end of the typical values observed in early studies but 
within the range found by Miettinen et al. (2011) in their IRDC cores. More evolved 
cores seem to show higher degree of ionization, which may be related to star-forming 
activities. 

4. In the four quiescent cores, the inferred characteristic scale for ambipolar diffusion is 
roughly 10~^ pc, and the couphng parameter of < 120 are within the regime of 
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marginal field-neutral coupling. The physical conditions may favor turbulent cooling 
flows. 

We thank Dr. H. Beuther and Dr. J. M. Rathborne for providing the 1.2 mm MAMBO 
and MAMB02 maps. We also thank Dr. H. Shang and Dr. R. Krasnopolsky for the helpful 
discussion on ambipolar diffusion. This research is supported by National Science Council of 
Taiwan through grants NSC 97-2112-M-007-006-MY3 and NSC 100-2112-M-007-004-MY2. 

A. Discussions on Individual Clouds 

G34.43+0.24.— The IRDC G34.43 {d = 3.7 kpc) contains nine cores (Rathborne et al. 
2006) with G34.43-MM2 being the most evolved core associated with the UC Hll region 
IRAS 18507+0121 of spectral type B0.5 (Mohnari et al. 1998; Shepherd, Niirnberger, & 
Bronfman 2004; Shepherd et al. 2007). The brightest millimeter core, G34.43-MM1, exhibits 
typical chemical signature of a hot molecular core (HMC) and has started internal heating 
with embedded source(s) (Rathborne et al. 2008). Additionally, G34.43-MM1, MM3, MM4, 
MM5, and MM8 are associated with extended Spitzer 4.5 /xm emissions, indicating possibly 
outflow activities (Chambers et al. 2009). 

IRAS 18151-1208.— The HMPO 118151 {d = 3.0 kpc) hosts four dusty cores (Beuther 
et al. 2002a) with MMl being the most evolved and dominant A'-band source, possibly 
driving II2 jets and CO outflows (Beuther et al. 2002b; Davis et al. 2004). By analyzing 
molecular line emissions and the spectral energy distributions (SEDs), Marseille et al. (2008) 
suggested an evolutionary sequence among three most compact cores: from the youngest 
I18151-MM3 perhaps in a pre-stellar phase, followed by MM2 as an HMSC with an embed- 
ded, mid-IR-quiet young protostar, to the most evolved MMl with mid-IR-bright protostars. 
Molecular outflows in I18151-MM1 and MM2 also suggest their more evolved stages (Mar- 
seille et al. 2008; Beuther & Sridharan 2007). 

IRAS 18223-1243.— The HMPO 118223 {d = 3.7 kpc) harbors a few dusty cores in 
a fllamentary structure with I18223-MM1 most evolved as a HMPO and others as HMSCs 
(Sridharan et al. 2005). Large SiO hnewidths suggest outflow activities in I18223-MM3 and 
MM4 (Beuther & Sridharan 2007). 
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Table 1. Properties of the Observed Cores 



Core 


R.A. 


Decl. 


a 


-^1.2 mm 


2R^ 


log 


Outflow^ 


Remark*^ 




(J2000) 


(J2000) 


(K) 


(Jy) 


(pc) 


(cin~^) 


(Y/N) 




G34.43+0.24 (IRDC; d = 3.7 kpc) 












MMl 


18 53 18.0 


+01 25 24 


18.5 


4.01 


0.19 


6.44 


Y 


HMC 


MM2 


18 53 18.6 


+01 24 40 


18.8 


4.33 


0.42 


5.42 


Y 


UC Hll 


MM3 


18 53 20.4 


+01 28 23 


15.5 


1.02 


0.38 


4.95 


Y 




MM4 


18 53 19.0 


+01 24 08 


17.6 


0.86 


0.38 


4.87 


Y 




MM5 


18 53 19.8 


+01 23 30 


14.3 


2.24 


0.89 


4.65 


Y 




MM6 


18 53 18.6 


+01 27 48 


14.0 


0.43 


0.62 


4.41 


N 




MM8 


18 53 16.4 


+01 26 20 


17.2 


0.36 


0.52 


4.44 


Y 




MM9 


18 53 18.4 


+01 28 14 


13.9 


0.53 


0.67 


4.41 


N 




IRAS 1^ 


U51-1208 (HMPO; d = 


3.0 kpc) 












MMl 


18 17 58.0 


-12 07 27 


20.8 


3.6 


0.25 


5.94 


Y 


HMPO 


MM2 


18 17 50.4 


-12 07 55 


17.8 


2.6 


0.34 


5.62 


Y 


HMSC 


MM3 


18 17 52.2 


-12 06 56 


16.0 


0.9 


0.40 


5.05 


N 




IRAS 18223-1243 (HMPO; d = 


3.7 kpc) 












MMl 


18 25 10.5 


-12 42 26 


17.5 


2.5 


0.55 


4.88 


Y 


HMPO 


MM2 


18 25 09.5 


-12 44 15 


15.1 


0.6 


0.47 


4.88 


N 


HMSC 


MM3 


18 25 08.3 


-12 45 28 


16.2 


0.8 


0.31 


5.42 


Y 


HMSC 


MM4 


18 25 07.2 


-12 47 54 


15.5 


0.3 


0.52 


4.43 


Y 


HMSC 



^ Cas temperature based on NH3 observations (Sakai et al. 2008). 

^Integrated 1.2 mm flux density and the deconvolved core size, deflned as the ge- 
ometric mean of the major and minor FWHMs. (Beuther et al. 2002a; Rathborne 
et al. 2006). 

*^Compiled from literature using Spitzer 4.5 fira emission and molecular outflows 
(Beuther et al. 2002b; Chambers et al. 2009; Lopez-Sepulcre et al. 2011; Marseille 
et al. 2008). 

'^Classiflcation from previous studies (Rathborne et al. 2006, 2008; Sridharan et al. 
2002, 2005). 
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Table 2. Parameters of Spectral Fits for N2H+ and N2D+ 



Core 


(TN2H+ 




A^fNoH+l 




Av 




A 


' max 




(mK) 


(mK) 


(10^2 cni-2) 


(km s~^) 


(km s""*") 








G34.43+0.24 
















MMl 


44 


9.5 


17 ±2 


57.550 ± 0.008 


4.81 ± 0.06 


(39 ± 6) X 10-^ 


3.7 


0.77 


MM2 


20 


7.7 


13 ± 1 


57.579 ± 0.005 


4.54 ±0.05 


(80 ± 7) X 10-^ 


9.8 


0.61 


MM3 


20 


7.8 


6.0 ±0.6 


59.66 ±0.01 


4.20 ±0.03 


0.007 ±0.002 


3.1 


0.32 


MM4 


21 


8.8 


5.8 ±0.5 


57.56 ±0.01 


4.76 ±0.03 


0.022 ±0.002 


4.2 


0.27 


MM5 


19 


7.1 


1.8 ±0.3 


57.89 ±0.03 


3.45 ±0.06 


0.033 ±0.005 


1.6 


0.12 


MM6 


23 


10.5 


0.59 ±0.09 


58.27 ±0.06 


2.3 ±0.1 


0.11 ±0.02 


0.8 


0.05 


MM8 


24 


8.5 


0.33 ±0.04 


57.2 ±0.1 


3.4 ±0.3 


0.00 ±0.02 


1.2 


0.02 


MM9 


17 


8.9 


1.6 ±0.2 


58.89 ± 0.02 


3.04 ± 0.06 


0.058 ± 0.007 


2.5 


0.11 


IRAS 1^ 


^151-120^ 
















MMl 


25 


10.6 


3.9 ±0.4 


33.203 ±0.008 


2.94 ±0.03 


0.010 ±0.002 


6.0 


0.27 


MM2 


22 


7.9 


6.6 ±0.8 


29.706 ±0.007 


4.01 ±0.04 


0.019 ±0.001 


3.9 


0.36 


MM3 


20 


8.3 


1.1 ±0.3 


30.67 ±0.02 


2.40 ±0.06 


0.064 ±0.007 


2.4 


0.10 


IRAS 1^ 


^223-1243 














MMl 


30 


8.6 


3.5 ±0.4 


45.41 ±0.02 


3.69 ±0.05 


0.021 ±0.002 


3.4 


0.20 


MM2 


38 


7.9 


1.9 ±0.3 


45.20 ±0.05 


3.7±0.1 


0.033 ±0.005 


1.1 


0.11 


MM3 


40 


8.6 


4.2 ±0.4 


45.54 ±0.03 


4.42 ±0.07 


0.021 ±0.002 


1.5 


0.21 


MM4 


43 


8.4 


0.9 ±0.2 


45.93 ±0.07 


2.5 ±0.2 


0.015 ±0.009 


1.2 


0.08 
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Table 3. Parameters of C^^O Spectra and Derived Quantities 



Core 


CTciSo 








>S'l.2mm 


fD 


logXe 




(mK) 


(Kkms-^) 


(10^^ cm-2) 


(K) 


(Jy/(34")') 






G34.43+0.24 














MMl 


51 


mi -1- n 1 


/I 71 -1- n no 




fi 17 
0.1/ 






MM2 


66 


23.8 ±0.1 


11.2 ±0.2 


> 34b 


4.97 


< 6± 1 


-5.30 


MM3 


63 


5.1 ±0.1 


2.32 ± 0.06 


32 


1.35 


9±2 


-5.27 


MM4 


77 


13.0 ±0.1 


6.0 ±0.1 


32 


2.74 


7± 1 


-5.77 


MM5 


57 


6.64 ±0.08 


3.0 ±0.1 




0.92 


14 ±3 


-6.01 


MM6 


52 


5.37 ±0.09 


2.44 ±0.09 




0.67 


13 ±3 


-6.62 


MM8 


55 


5.2 ±0.1 


2.42 ±0.07 




0.72 


10 ±2 




MM9 


54 


4.68 ±0.08 


2.12 ±0.07 




0.75 


16 ±3 


-6.28 


IRAS 1^ 


U51-120; 


8 












MMl 


54 


16.13 ±0.09 


7.8 ±0.2 


27.3 


2.41 


6± 1 


-5.37 


MM2 


59 


12.0 ±0.1 


5.6 ±0.1 


19.4 


1.85 


10 ±2 


-5.71 


MM3 


49 


8.33 ± 0.08 


3.8 ±0.2 




0.81 


8±2 


-6.33 


IRAS 18223-1243 












MMl 


69 


20.3 ±0.1 


9.4 ±0.2 


31 


1.74 


3.7 ±0.7 


-5.80 


MM2 


60 


10.10 ±0.09 


4.6 ±0.1 




0.54 


6± 1 


-6.03 


MM3 


72 


8.9 ±0.1 


4.1 ±0.1 


18 


0.71 


7± 1 


-5.77 


MM4 


73 


4.6 ±0.1 


2.09 ±0.09 




0.24 


6± 1 


-5.65 



^ Dust tcperature based on SED fits in literature (Rathborne et al. 2005; Marseille 
et al. 2008; Bcuther et al. 2010). 

^ Dust temperature adopted from G34.43-MM1. 
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Fig. 1.— (a)-(h) Spectra of N2H+ (3 - 2) shifted by 8 K (top histogram) in the order 
of G34.43-MM1, MM2, MM3, MM4, MM5, MM6, MM8, and MM9, superposed with the 
thirty-eight hyperfine component synthesis spectra (sohd curves) assuming a single excita- 
tion temperature based on NH3 observations (Sakai et al. 2008). The frequency of each 
individual hyperfine component in the model is labelled with a short bar on the top. Spectra 
of N2D+ (3 — 2) multiplied by 10 and shifted by 6 K (middle histogram) toward the corre- 
sponding cores superposed with synthesis spectra (solid curves) and spectra of C^^O (2 — 1) 
(bottom histogram) are also shown. 
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Fig. 2.— (a)-(c) Similar to Fig. [T] but for I18151-MM1, MM2, and MM3. Spectra of 
N2H+ (3 — 2) are shifted by 9 K and those of N2D+ (3 — 2) are multiphed by 10 and shifted 
by 7K. 
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Fig. 3.— (a)-(d) Similar to Fig.|T|but for I18223-MM1, MM2, MM3, and MM4. Spectra of 
N2H+ (3 - 2) are shifted by 10 K and those of N2D+ (3 - 2) are multiphed by 10 and shifted 
by 8 K. 
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Fig. 4. — (a) Deuterium fractionation, Dfrac, vs. gas temperature, Tg, shows a general 
decreasing trend. This correlation has Spearman's p rank correlation coefficient of p = —0.67 
with significance of 0.6 % and Kendall's r rank correlation coefficient of r = —0.50 with 
significance of 1.0 %. Cores in G34.43, 118151, and 118223 are marked by ffiled circles, open 
diamonds, and crosses, respectively. Colors differentiate objects with known categories: blue 
for one UC Hll region, green for one HMC and HMPOs, and red for HMSCs; black color 
reserves for objects that were not previously classified, (b) -Dfrac vs. fitted linewidth, Av, 
exhibits a decreasing trend that has p = —0.61 with significance of 1.6 % and r = —0.49 
with significance of 1.2%. (c) -Dfrac, vs. CO depletion factor, fi). Excluding G34.43-MM1 
that has an unusually large fn, the distribution of all other cores shows a weaker increasing 
trend that has p = 0.49 with significance of 7.8 % and r = 0.34 with significance of 9.0 %. 
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Fig. 5. — (a) Linewidth N2H+ (1 — 0) spectra from Sakai et al. (2008), Av{X), vs. linewidth 
of our N2H+ (3 — 2) spectra, Av. The angular resolutions for the J = 1 — and 3 — 2 
observations are 18" and 27", respectively. All the cores have larger linewidth in the J = 3 — 2 
transition. Cores in G34.43, 118151, and 118223 are marked by filled circles, open diamonds, 
and crosses, respectively. A solid hne gives the locus of equal linewidths between the two 
spectra. The ordinate label, Av{X), stands for the linewidth of the line shown on the 
top. (b-d) Linewidth of NH3 (1,1), (2,2), and (3,3) spectra from Sakai et al. (2008) vs. 
hnewidth of our N2H+ (3 — 2) spectra. The angular resolution for NH3 observations are 
roughly 73". Except in the most quiescent cores, G34.43-MM6 and I18151-MM3, linewidths 
of the N2H+ (3 — 2) spectra are larger than those of NH3 (1, 1) and (2,2) spectra. On the 
other hand, the NH3 (3, 3) transition with higher E^p = 124.5 K seems to tracer warmer 
component with linewidth larger than N2H"'" (3 — 2). Note that only ten cores, not including 
any of the four quiescent ones, were detected in NH3 (3, 3) emission (Sakai et al. 2008). 



